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Abstract 

We study planar nematic equilibria on a two-dimensional annulus with strong and weak 
tangent anchoring, within the Oseen-Frank and Landau-de Gennes theories for nematic liquid 
crystals. We analyse the defect-free state in the Oseen-Frank framework and obtain analytic 
stability criteria in terms of the elastic anisotropy, annular aspect ratio and anchoring strength. 
We consider radial and azimuthal perturbations of the defect-free state separately, which yields a 
complete stability diagram for the defect-free state. We construct nematic equilibria with an arbi¬ 
trary number of defects on a two-dimensional annulus with strong tangent anchoring and compute 
their energies; these equilibria are generalizations of the diagonal and rotated states observed in 
a square. This gives novel insights into the correlation between preferred numbers of defects, 
their locations and the geometry. In the Landau-de Gennes framework, we adapt Mironescu’s 
powerful stability result in the Ginzburg-Landau framework (P. Mironescu, On the stability of 
radial solutions of the Ginzburg-Landau equation, 1995) to compute quantitative criteria for the 
local stability of the defect-free state in terms of the temperature and geometry. 


1 Introduction 

Nematic liquid crystals (LCs) are classic examples of partially ordered materials that combine the 
fluidity of liquids with the orientational order of crystalline solids mm- Nematics have generated 
substantial scientific interest in recent years because of their unique optical, mechanical and rhe¬ 
ological properties [3] and notably, nematics form the backbone of the multi-billion dollar liquid 
crystal display (LCD) industry. Defects are a key feature of nematic spatio-temporal patterns in 
confined geometries. A lot remains to be understood about the structure of defects and how they 
can be created, controlled and manipulated to yield desired properties. In this paper, we revisit 
the classical problem of a nematic sample in a two-dimensional (2D) annulus with strong or weak 
tangent boundary conditions separately and no external fields. Our work is motivated in part by re¬ 
cent experiments mis], on rod-like /d-virus particles within shallow, annular microscopic chambers. 
Within these chambers, multiple states are observed, including a radially invariant defect-free state 
and states with regularly arranged defects on the boundary, see Figure 

The model problem of a defect-free state in annular wells has received a great deal of attention 
in the past, especially within the Oseen-Frank (OF) theory for liquid crystals. Here, we give a brief 
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overview. The papers la El El El uni m [12] are directly relevant to our work. In El) th® author 
studies the stability and multiplicity of nematic radial equilibria on a 2D annulus with strong uniform 
anchoring. In [HlEllII]) tbe authors approach the same problem with a more applications-oriented 
perspective motivated by the classical Freederickzs transition. They study confined nematic samples 
between two concentric cylinders with weak anchoring on one lateral surface and strong anchoring on 
another, subject to an external magnetic field. The authors primarily consider the stability of three 
characteristic conhgurations, referred to as ‘radial’, ‘azimuthal’ and ‘uniform’ and obtain explicit 
estimates for the critical threshold held in terms of the OF elastic anisotropy, cylindrical aspect 
ratio and anchoring strength. The defect-free state in our paper is analogous to the azimuthal state 
on an annulus, with tangent boundary conditions on both circular boundaries [3 El uni E] • Whilst 
our results in the OF case, both with strong and weak anchoring, are partially captured by the 
results in ElElElEnilll], our method of proof is different. We compute the second variation of the 
anisotropic OF energy and study the resulting eigenvalue problem directly. We compute stability 
curves in terms of the elastic anisotropy 6, the annular aspect ratio b, the anchoring strength, a, 
on both boundaries and the order of the azimuthal perturbations, k. Such stability diagrams can 
provide useful insight into stabilization and de-stabilization effects in simple geometries, i.e. they 
can quantify the response of the defect-free state to different types of azimuthal perturbations. For 
k = 0, our results reduce to the previously reported results in El SISl EHl III]- Finally, in [6], the 
authors demonstrate that the planar radial state loses stability, when the ratio of the inner radius to 
the outer radius is smaller than a critical value, and the planar radial state escapes into the ‘third’ 
dimension. In the authors study nematic samples confined between two co-axial cylinders, with 
emphasis on higher-dimensional biaxial effects, which are outside the scope of the present paper. 

In ESI E] , the authors study nematic equilibria conhned to shallow square and rectangular wells, 
subject to strong tangent boundary conditions on the edges. They construct analytic approximations 
for the experimentally observed ‘diagonal’ and ‘rotated’ states and obtain explicit expressions for 
their respective energies. The diagonal state always has lower energy than the rotated states. We 
take this work further by constructing explicit solutions of the Laplace’s equation on an annular 
sector of angle 2iTjN (with iV G N) with tangent boundary conditions. These ‘sector’ solutions 
mimic generalized ‘rotated’ and ‘diagonal’ nematic equilibria on 2D annuli, featuring a total of 2N 
defects pinned to the inner and outer boundaries. We compute the respective energies in the OF 
framework and hnd that there are interesting energetic cross-overs, in terms of the number of defects. 
Within a sector, we find that the rotated states typically have lower energies than the corresponding 
diagonal states for small values of A^. As N increases, the sector approaches a rectangle and the 
generalized diagonal state becomes energetically preferable, consistent with the energetic trends in 
a rectangle. States with defects on the boundary of an annular well are considered in |15] . The 
authors work within the one constant Oseen-Frank framework with weak anchoring and use an 
expression for the director based on a linear combination of defect-type solutions located at equally 
spaced locations on the boundary. This expression, although not a solution of the Euler-Lagrange 
equations, allows the authors to identify possible director profiles and obtain estimates to the surface 
and elastic energies. These energies demonstrate that states with boundary defects may have lower 
energy than the defect-free state in certain parameter regimes by including weak anchoring. In this 
paper, we demonstrate that in the strong anchoring regime, this is only possible for liquid crystals 
with large elastic anisotropy. 

Equilibria with boundary defects can be interesting in physical situations, particularly with 
Neumann boundary conditions or when we simply specify topological degrees on the boundaries, as 
opposed to Dirichlet conditions. For example, in m , the authors study 2D vector fields, u = {ui,U 2 ), 
on a multiply-connected 2D domain with prescribed topological degrees on the outer boundary and 
on the inner boundaries enclosing the ‘holes’. The authors study minimizers of the Ginzburg- 
Landau functional on such domains and hnd that the inhmum energy is not attained. Minimizing 
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sequences develop vortices or boundary defects in certain asymptotic limits. Our analysis of OF 
equilibria with boundary defects could be instructive for such model LC problems, with topological 
degree boundary conditions; in fact, they give quantitative information about minimizing sequences. 
Equally importantly, our analysis is directly relevant to recent experimental work on colloidal samples 
in shallow annular wells wherein the authors observe states with defects pinned to the lateral surfaces 
015], such as those shown in Figure[^ The relative energies of the generalized ‘diagonal’ and ‘rotated’ 
states give qualitative insight into the relative observational frequencies of the experimental states 
that exhibit boundary defects. 



Figure 1: Nematic director fields identified experimentally in 015]. 

We complement our analysis in the OF case with work within the more general and powerful 
Landau-de Gennes (LdG) theory for nematic liquid crystals. In [1^, the author studies the generic 
+1-degree vortex within the Ginzburg-Landau (GL) theory for superconductivity and derives a 
powerful stability result for the -|-1-degree vortex on a 2D disk (that contains the origin), with 
Dirichlet radial boundary conditions. In 2D, the LdG theory for nematic liquid crystals reduces 
to the GL theory and the nematic state is fully described by a scalar order parameter and a two- 
dimensional vector held, referred to as ‘director’ in the literature m- We dehne the nematic 
defect-free state in the LdG framework, in terms of a minimizer of an appropriately dehned energy 
functional, and adapt Mironescu’s proof to a 2D annulus. We address technical differences due to 
the change in geometry and demonstrate local stability in the low temperature limit, and derive a 
stability criteria in terms of temperature and geometry. 

The paper is organized as follows. In Section]^ we review the OF and LdG theory for nematic 
LCs. In Section and we focus on the OF theory and compute stability criteria for the defect-free 
state as a function of 6, b and a as introduced above. In Section]^ we construct generalized ‘diagonal’ 
and ‘rotated’ states in an annulus with an arbitrary number of boundary defects and compute the 
corresponding energies. Finally, in Section we compute criteria for the local minimality of the 
LdG defect-free state on a 2D annulus as described above and conclude in Section 0 with future 
perspectives. 

2 Theory and Modelling 

Nematic liquid crystals are complex liquids wherein the constituent rod-like molecules move freely 
as in a conventional liquid and tend to align along certain preferred directions 0(121 US). The 
Landau-de Gennes (LdG) theory is one of the most general continuum theories for nematic LGs to 
date mm- The LdG theory describes the nematic state by a macroscopic order parameter, known 
as the LdG Q-tensor, that is a macroscopic measure of the LG anisotropy. For three-dimensional 
problems, the LdG Q-tensor is a symmetric, traceless 3x3 matrix with five degrees of freedom 0 [T9| . 
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A nematic phase is said to be (i) isotropic if Q = 0, (ii) uniaxial if Q has two degenerate non-zero 
eigenvalues and (iii) biaxial if Q has three distinct eigenvalues. The LdG theory is a variational 
theory and a prototypical LdG energy functional is of the form 

-^[Q] := wiQ, VQ) -b fsiQ) dn where /b(Q) = J (1) 

w (Q,VQ) is an elastic energy density that penalizes spatial inhomogeneities and fs is a bulk 
potential that drives the nematic-isotropic phase transition as a function of the temperature or 
concentration. In Equation Q, A < 0 and is proportional to the rescaled temperature and B,C > 0 
are positive material-dependent constants. Whilst studying 2D nematic equilibria on 2D domains, 
the LdG Q-tensor reduces to a symmetric, traceless 2x2 matrix of the form 


Q 


s{x,y) 



( 2 ) 


where n is a two-dimensional unit-vector and I 2 is the identity matrix in 2D [TS]. It is clear from 
Q that a 2D Q-tensor only has two degrees of freedom: the real-valued scalar order parameter, 
s{x,y), that is a measure of the degree of the orientational order and the unit-vector field, n G S\ 
that represents the distinguished direction of alignment of the nematic molecules. For Q as in Q, 
trQ^ = 0 and hence, the LdG energy reduces to 


I[Q] := ^ w{Q, VQ) + ^trQ2 + ^ (trQ^)' dD. 


(3) 


We work with low temperatures, described by A < 0, and the one-constant elastic energy density, 

, so that the LdG energy in (3) reduces 


Qq 


t(;(Q, VQ) = ^|VQf. We introduce the scaling Q = 
to the Ginzburg-Landau functional in superconductivity 


/[Q]= [ |VQ|2 + ^(|Q|2-i)2da (4) 

J O 

As in any problem in the calculus of variations, the problem of studying nematic equilibria is math¬ 
ematically equivalent to a study of local and global energy minimizers. In particular, we obtain an 
explicit relation between geometry and temperature that guarantees local stability of the defect-free 
state. 

The Oseen-Frank (OF) theory is a simpler continuum theory for nematic liquid crystals restricted 
to uniaxial nematic phases with a constant scalar order parameter m- In this case, the macroscopic 
order parameter is simply a unit-vector field, n, that defines the unique direction of molecular 
alignment. For a 2D problem, the OF energy functional is given by 


E[nl := ^ (n X (V X n))" dil, (5) 

where D C is a 2D domain and the elastic constants Ki and are associated with splay and bend 
director deformations [T]. For 2D domains, n can be conveniently written as n = (cos 0, sin 0,0), 
where 0 : —>• M is a function of the planar polar coordinates, (r, (j)), and the OF energy ([^ reduces 
to 

■= jj (^os{6 — 4>)—— s\n{6 — 4>)6^ -\—— (f))—+ cos{6 — dD. (6) 
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We define 6 := 1 — ^ to he the measure of elastic anisotropy. The corresponding Euler-Lagrange 
equations are: 


1 - 


+ (5 ( - sin(20 — 2(f)) 


‘^Sr(j> ^4) 

- + 72 “ “ 


+ ^ cos(20 - 2(1)) (orr - ^ ~ + 



= 0 , 


(7) 


subject to appropriate boundary conditions for 9 on dCl. 

In Sections 1^ and 1^ we focus on the defect-free state on a 2D annulus with tangent boundary 
conditions. The rescaled 2D annulus is defined by 


D = |(r, 6<r<l, t) < (f) < 2 t('^ , 

where the rescaled radius b := is the ratio of the inner and outer radii. 

-rtouter 

consider the reduced domain 


( 8 ) 


In Section 
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we 


Qn 


I (r, (f>) 


2 . 


b <r <1, 


0 < (/)< 



(9) 


where N € N, which allows us to consider more complex equilibria with defects. The tangent 
boundary conditions require that 


TT 


9(1, (f)) = 9{b,(f)) = 2 


( 10 ) 


Equation (10) is the Dirichlet version or the strong version of the tangent boundary conditions. The 
tangent boundary conditions can be weakly implemented too, as will be illustrated in Section]^ The 
OE defect-free state is simply defined by 


TT 


9*{r,(f)) = </>+ 2 ’ 


{r,(f)) G n, 


( 11 ) 


with no radial variations. A simple computation shows that E[9*] = vrifTslog (^) , so that the OF 
energy diverges in the limit 6—7-0. In the LdG framework, by analogy with similar work on degree 
-|-1-vortices in the GL theory |20] . we define the LdG defect-free state to be 


Q*{r,(j)) 



( 12 ) 


where 

n* = (cos0*,sin0*,O) 

and s is the unknown scalar order parameter, independent of cj). In what follows, we study the 
defect-free state in both theoretical frameworks and make comparisons to competing states with 
boundary defects. 


3 The OF Defect-Free State with Strong Anchoring 

We firstly consider the OF defect-free state in 0 with Dirichlet boundary conditions. Our result 
is similar to the stability estimate in mm but our method of proof is different. We consider 
the second variation of the OF energy about 9*. The strict positivity of the second variation is a 
sufficient criterion for local stability, whereas 9* is unstable if the second variation is negative for 
some admissible perturbation |21j . 
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We perturb 0* by er/, where r/( 6 , 4>) = 77 ( 1 , (^) = 0 in accordance with the Dirichlet conditions in 
(10). A straightforward computation shows that the second variation is given by 


,2^__d^E[e* + er^\ 


:= 


de'^ 


e=0 


= JJ^{Ki-Ks){^ + Vry+ K3{Vrjf dn 
= K3 ff {Vvf -s{y + vry dn 


(13) 

(14) 


and it follows immediately that 6‘^E > 0 for S < 0 and any non-trivial rj. 

Proposition 1: For any admissible 'r]{r), there exists drj < 1 such that S'^E < 0 for <5^ < <5 < 1. 
Proof. We define the functional E 2 , 

2 


E 2 iv, 6 ,b):= - 6 (y-+rjr) dA. 

For rj independent of E 2 simplifies to 

E 2 {ri{r), 6 ,b) = 2 tt | (r?r)^ - ^ rdr. 


(15) 


(16) 


For a fixed b and E 2 is a continuous function of 5 which is bounded below by E 2 {r]{r), 1,6). 
Using the boundary conditions 7 /( 1 ) = 77 ( 6 ) = 0, this lower bound simplihes to 


E 2 {r]{r),l,b) = —27r [ — + 2 r]r]rdr =—271 [ —dr 

Jb r Jf, r 


(17) 


and hence E 2 {r]{r), 1,6) <0 for all rj and any 0 < 6 < 1. We use the continuity of E 2 with respect 
to 6 to conclude that for all 77 (r), there exists 6 ,; < 1 such that < 0 for all 6 ,^ < 6 < 1 . □ 

The preceding observation demonstrates that the defect-free state loses stability on a 2D annulus, 
for certain choices of ( 6 , 6 ). We make this observation more rigorous by studying the Sturm-Liouville 
problem associated with the minimization of 6'^E as shown below. The corresponding Euler-Lagrange 
equation is: 

dr' 


bd I d(t>lt> I /I e'l / \ \ _ O 


(18) 


where 77 = 0 on r = 6 ,1 and 77 is 27r-periodic. Without loss of generality, we can consider separable 
solutions of the form rj = e*^‘^/fc(r), where /c G Z. The function fk{r) is a solution of 


r/^(r) +r 2 /"(r) + ^ 3 ^/fc(?^) = 0 , 


(19) 


with boundary conditions /fc( 6 ) = /fc(l) = 0. Equation (19) is an Euler ordinary differential equation 

( 20 ) 


and we take fk to be of the form fk = where m and k are related by 

k ‘^-6 


mf = 


1-6 


For k = 0, the general solution of ( 19|) is given by 
/o(r) = A sin 


1-6 


log(r) -|- B cos 


1-6 


log(r) 


( 21 ) 
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and applying the boundary conditions /o(l) = fo{b) = 0, we find that 


foA'^) = ^sin 7 rn 


loArA 

iog( 6 )y 


where Sn = 


vr^n^ 


TT^n^ _|_ iog(5)2 


( 22 ) 


For A; > 1, the boundary conditions on naturally rule out any non-trivial solutions. The case A; = 0 
is a typical Sturm-Liouville problem [22] with eigenfunctions {/o,n} and corresponding eigenvalues, 

{Jn}- Similarly, substituting rj = A sin ^tt^ into (13), we find that the corresponding second 

variation is negative, 6 '^E < 0 for <5 > (ii. Therefore, the defect-free state loses stability along the 
curve, <5 = (5i (6) defined by 


<^1 = 


TT 


TT^ -|- log( 6)2 ' 


(23) 


This stability criterion (23) is known in the literature mra]. 

We next establish that the defect-free state undergoes a supercritical pitchfork bifurcation at 
above. We let -|- €^62 and consider perturbations of the form 


9 = 6* + erji + e^??2 + + • • • 


(24) 


where the functions r]i vanish on r = 6,1 and are 27r-periodic in cj). We substitute (24) into the 


Euler-Lagrange equation ([^ and compare terms of order e. We define the operator L to be 




Cn :=AA + J^ + ^l_Si){Vrr + ^). 


Vr 


(25) 


A straightforward computation shows that Crji = Crj 2 = 0. Recalling (18), we deduce that rji and 
772 are proportional to sin ( io|(b) ) • Comparing terms of order e^, we find 


where the function R is defined to be 


£773 = R(r), 


(26) 


R{r) AhM') 

2 r^\T?{h) V la( 6 ) y 


+ lA{b)) 

6 r^ ln^(6) 




(27) 


From the Fredholm Alternative Theorem, the differential equation (26) admits a solution if and only 
if the following compatibility condition holds 


sm ) R{r)rdr = 0, 


This condition reduces to 


log(6) 


AiAHi - 262) = 0 , 


(28) 


(29) 


and therefore the amplitude. A, is only defined for <^2 > 0 and we have a supercritical pitchfork 
bifurcation at <5 = (5i. 

Finally, we explicitly construct a solution of the Euler-Lagrange equations 0 that gives us 
insight into how the defect-free solution deforms for <5 > and the structure of the corresponding 
global energy minimizers for <5 > 5i. We seek an exact solution of the form 


e = 9* + U{5,r), 


(30) 
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where U satisfies C/(5,1) = U{6,b) = 0. We consider functions with at most one critical point, 
solutions with more critical points exist but have a higher energy. We dehne t = — log(r), substitute 
(30)into 0 and then U satisfies 


((5 cos^{U{t)) — 1) 






+ 1 = 0 , 


(31) 


subject to U{6,0) = 0 and U{6,log{l/b)) = 0. 
We can define U{S, t) implicitly by 


/ 1 — 5 cos^ u 

Iq y 5 cos^ u — 5 cos^ Uq 


du = t, 


(32) 


where Uq is the maximum value of U{5,t). 

In Figure]^ we plot the spiral-like solution (30) for b = 0.2 and 5 = 0.95. Numerical computations 
of the energy show that this spiral-like solution has lower energy than the defect-free solution in the 
parameter regime > 5i. 
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Figure 2: The director field for 9 = 9* + U, with b = 0.2 and 5 = 0.95. 


When (5 = 1, we can solve (31) exactly to find 

= arccos -- 

^ ^ \b + l 


exp(i) + exp(-t) 


(33) 


with associated energy, E = 2tt • Given that (33) is an exact solution, we can study its stability 

in terms of the second variation; 


:= + (1 - g‘^)rif if 


2g‘^--ig^+g'^ + l + 2gg" -2g-^g‘ 


3 ..'/ 


1-g^ 


dll. 


(34) 


where g{t) = exp(t) -|- exp(—t). The function 5 < 1 for t G (0, — log( 6 )) and the coefficient of 

simplifies to 

26 - 62 _ 1 ^ 

62e2t + e-2t _ 52 _ 1 ^ 1- ( 35 ) 

Therefore, the second variation is positive for all admissible perturbations and we have established 
local stability of the spiral-like solution ( |30[ ) at 6 = 1. In summary, we have constructed an exact 
solution of the Euler-Lagrange equations, (1^, which is a radial variant of the defect-free state and 
is locally stable for 5 close to unity. 





















4 The OF Defect-Free State with Weak Anchoring 


We add a surface energy to the Oseen-Frank energy in ([^ , that enforces preferential planar anchoring. 
Our approach is similar to the work in m and our main contribution is the stability diagram in 
Figure that is particularly useful in quantifying the interplay between anisotropy, geometry and 
surface effects in the stability of the defect-free state. We employ the Rapini-Papoular surface energy 



where 90 consists of two concentric circles with radii r = b and r = 1, kF is an anchoring coefficient 
and ds is an arc-length element along 90. We define a := to be the ratio of the 

outer radius to the extrapolation length d, and the three key parameters are 6, b and a. The 
corresponding Euler-Lagrange equation is 0 with boundary conditions 


2-5 

2 


-k 


<5 ( 0 , 


j 

— sin(20 — 2(j)) + Or cos(20 — 20) 


(y 

+ — sin(20 — 20o) = 0, 


(37) 


on r = 1 and 


2-5 


sin(26' - 20) -k Or cos(26l - 20) j -k | sm(20 - 20o) = 0, (38) 


on r = 5. It is straightforward to check that the OF defect-free state, 0*, defined in ( |11| ), is a solution 
of 0 with the above boundary conditions. By analogy with our work in Section]^ we compute the 
second variation as shown below: 


5^i;[0*] := jj {Vr]f - 5 {^ + + a j r] 


‘ds 


(39) 


where r/ is a perturbation about 6*. Without loss of generality, we assume rj = /fc(?’)e*^'^, where 

we refer to k as being the azimuthal order of the perturbation. For each k, the optimal fk is a 
solution of 

rfkir) + r'^fkir) + — ^fk{r) = 0 , ( 40 ) 


with boundary conditions 


l-<5 


fkir) = on r = 1, 

fkin = 


1-5 


fk{r) on r = b. 


( 41 ) 

( 42 ) 


Standard computations, following a change of variable t = log (7), show that fk{t) is given by 



for A: = 0 


for A: > 1 


( 43 ) 

( 44 ) 


and the boundary conditions require the following ‘compatibility condition’ between a, 5, b and k: 
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tan 



ail + b)^6{l-6) 
ad — ab6 + — 6 


for k = 0 


(45) 


tanh 



(1 — 6)a{l + b) 


5a + a^b — ab5 — 5 + /c^(l — 6)\ 1 — 5 


k^-5 


= 0 


for A: > 1. 


(46) 


For a given k G N, there will typically be multiple solutions {5i^k, 52,k, ■ ■ ■} of the compatibility 
relation above. We are interested in the smallest solutions, whilst defining stability curves in 
the {5, Q:)-plane for different values of k and b. 


4.1 The case k = D 

The case A: = 0 is contained in the work of m- However, our method of proof is different and of 
independent interest. We first point out that for k = 0, the boundary-value problem (40) - (42) 
is a Sturm-Liouville problem with a complete set of eigenfunctions and corresponding eigenvalues, 
{<^ 1 , 0 ) ^ 2 , 0 , • ■ ■) ^n,o} l22]. The eigenvalues, are solutions of the compatibility condition 

ail + b)y^5{l-5) 


tan 



a5 — ab5 + a^b — 5 


= 0 , 


(47) 


which has an infinite number of solutions for any fixed triplet (a, 5, b). The corresponding eigenfunc¬ 
tions are given by 


flnflir, <t>) = ^sin 



+ 41 


\/ <5n,o(l <^n,o) 


a — 5, 


COS 


n,0 



(48) 


In Figure]^ we compute as a function of b, for different values of a. As a —>■ oo, the compatibility 
condition becomes tan ^y^^^^ln (^)^ —)• 0 and strong anchoring 


result. 



Figure 3: The first eigenvalue for fixed a and the strong anchoring limit. 


We further investigate the loss of stability of 0* at 5 = 5i^o- Let 0 = 0*+ erji + e^r ]2 + and 
^ = <^ 1,0 + as before. Then r]i and rj 2 are p roportional to Equation (48) with n = 1. The 
corresponding compatibility condition for 773 to exist is 


mhR{r)rdr = rrj^h 




— + v[ 

r 


+ 3^ I « - TT 


rr]3h (^<52 + r/'i) - + 


d+o 

r 

<^ 1,0 

~V 


c 2 / 

- difiViVi 


r=l 


(49) 
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where rj^h is the solution of the homogeneous problem as given below: 


V 3 hir,(p) = sin 




, V <^1,0(1 — < 5 i,o) 

+ '■ .. cos 

a — Oi^o 



and 


2 // , 


R{r) ^62\^- + Vi- - <^1,0 [viVi - ^ + + 


The compatibility condition yields 


A^E 3 - 52AE1 = 0, 


( 50 ) 


(51) 


(52) 


where £"3 and £1 are functions of ( 5 i^ 0 ) b and a. We numerically can evaluate £1 and £3 at 
^ = ^ 1 , 0 ( 0 ^)^) to find that £1 and £3 are both positive. Hence, we have a supercritical pitchfork 
bifurcation at J = 61 ^ in the weak anchoring framework. 

4.2 The cases k > 1 


For k = 1, the compatibility condition (46) simplifies to 

(1 + b){lP‘a^ — ab^ + 26ab — 26b + b — a^b — a + 26 — 1 ) 


with a unique solution 


( 6 ^ + l)(5a + a^b — Sab — 25 + 1 


1 ab'^ + a^b + a + 1 — b'^a'^ — b 

= 2--■ 


= 0 , 


(53) 


(54) 


We note that G (0,1) iff 0 < a < 1 and that (5iq > 0.5 for all b. 


For k > 1, the compatibility condition (46) has solutions if and only if 


(1 — (5)a;(l + b) 


k ‘^-6 


6 a + a^b — abS — (i + ^^(1 — 6 ) \ 1 — 5 


In particular, this implies < 6 i^k < 1) which in turn requires that a < 1. Therefore, 

the compatibility condition (46) only admits azimuthal perturbations with A: > 1 if 0 < a < 1. In 


< 0 . 


(55) 


Figure we plot the stability curves in the (5, Q:)-plane, for hxed values of k and b. These curves 
define the stability of the defect-free state with respect to the different kinds of perturbations e.g. 
the case A: = 0 corresponds to purely radial perturbations whereas A: > 1 describe perturbations with 
azimuthal dependence. As 6 increases, the system loses stability upon crossing one of the curves 
and as is evident from Figure larger values of a stabilize the defect-free state with respect to 
all perturbations. Similar comments apply to the confinement parameter b. As b increases, the 
defect-free state is stable with respect to larger class of perturbations, as can be seen from Figure]^ 
The critical value a = 1 corresponds to ^ = £outer and the defect-free state is stable with respect to 
all azimuthal or symmetry-breaking perturbations for a > 1 or ^ > £outer- 


5 Nematic Equilibria with Defects 


In this section, we consider nematic equilibria on a 2D annulus with regularly spaced boundary 
defects in the Oseen-Frank framework. 

To model states with regularly spaced defects along cAI, we split D into N sectors with £ G N, 
and define the region to be 


Q.N 


I (r, (/)) G : 6 < r < 1, 0 < (?!) < 



(56) 
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(a) b = 0.1 (b) b = 0.5 



Figure 4: The smallest solution to the compatibility condition (45) and the strong anchoring limit, 
and the only solution to (46) for A: = 1, 2, 3 in the ((5, a)-plane for fixed b. 


We assume strong tangent boundary conditions on all four edges of fljv, r- = 6 and r = 1, (/> = 0 
and (f) = ^, so that there are necessarily discontinuities at the four corners. We only consider states 
with boundary defects of strength m = ± 1 , states with defects of higher strength are expected to 
have higher energy [TJ.The m = +1 defect has a local radial splay profile and the m = — 1 defect 
has a local bend profile. There are 4 distinct arrangements of defects on the boundary, 912jv, as 
shown in Figure to be considered. Three of these arrangements correspond to generalizations of 
the rotated state found in a square and rectangular wells, denoted Ui, U 2 and U 3 , where the director 
field connects two adjacent + 1 -defects along an edge, and the fourth arrangement is a generalized 
diagonal state, D, for which the director field connects two diagonally opposite +l-defects. In the 
Ui state, the two +l-defects lie along the edge r = b; the two -|-l-defects lie at the vertices of the 
edge r = 1 for the rotated state U 2 and the two -|- 1 -defects are located at the vertices of the edge, 
^ = 0, for C/ 3 . Other rotated and diagonal states are rotationally equivalent to one of the four cases 
enumerated in Figure The overall configuration in 12 is a superposition of the different states 
in the different sectors, We note that for odd N, a superposition of D or C /3 states would be 
physically unrealistic as it would require the superposition of a -|-1 and a — 1 -type defect at one of 
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the vertices. 



(a) Rotated state 1 {Ui) (b) Diagonal state 1 {D) (c) Rotated state 2 {U 2 ) (d) Rotated state 3 (C/ 3 ) 


Figure 5: The four distinct defect arrangements with the sector r 24 , with the strength of the defects 
denoted by the numbers in the corners, and sketches of the director. 


5.1 The One Constant Approximation 

Using the one-constant approximation, for which Ki = K 3 = K, the OF energy reduces to 

m--= JJ^j\ve\‘^dA. (57) 

The key step is to construct solutions of the Laplace equation, = 0 on Oat, with Dirichlet 
boundary conditions (determined by the tangent conditions) on the four edges. The tangent condi¬ 
tions require that 0 = (^±|onr = 6 and r = 1; 9 = 0 on 4> = 0 and 9 = { ^ ± vr} on the edge 

(j) = Any admissible 0, subject to these Dirichlet boundary conditions, can be written as 

9 = aocf) + aifi -b 02/2 + 03/3 -b (58) 

where the canonical functions fi are solutions of the Laplace equation with boundary conditions: 

• /i(uO) = /i(r,27r/iV) = 0, /i(6,(/)) = 0 and /i(!,</>) = 1 

• / 2 (r, 0 ) = /2(r,27r/iV) = 0, f 2 {b,c/)) = 0 and / 2 ( 1 , 0 ) = </> 

• h{r, 0) = h{r, 2 'k/N) = 0, h{h, cj)) = l and /3(1, (/>) = 0 

• U{r, 0) = fi{r, 2 tt/N) = 0, fi{b, (j)) = (f) and /4(1, </>) = 0 . 

We compute the canonical functions, /i... / 4 , using separation of variables, to be: 


00 4 sin 
/i(u<(') = ^- 

n=l 


{2n-l)N 

2 


(2n-l)JV .. (2n-l)JV 

- 2 - —b-(2n-l)N^ - 2 - 

l_6-(2n-l)iV 


(2n — l)7r 


/2(U </') = X] 


fsir, 4>) = 


hirA) = 

n=l 


sin 


nN ,T nN 


n=l 


00 4sin(M^ 

V 

'^)( 

2 ^ 

n=l 

(2n 

00 4(-l)"+^sin( 

Nniv ^ 


Nn 


{2n-l)N {2n-l)N 

2 —r 2 

(2n-l)Af (2n-l)iV 

2 -b 2 


nN nN 

■“ 2 “— 
nN nN 


Nn 


(59) 


(60) 


(61) 


(62) 
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(a) Rotated state 1 (C/i) (b) Diagonal state 1 {D) 




(c) Rotated state 2 {U 2 ) 


(d) Rotated state 3 (t/ 3 ) 


Figure 6 : The director of the four states within Q 4 and b = 0.5, in the one constant approximation. 


Sample director plots using the canonical functions are shown in Figure for N = A and b = 0.5. 
Due to the discontinuities in 0 at the corners of Dtv, we regularize the domain by removing disks of 
radius e about each of the defects (to prevent the energy (57) from diverging) and the new regularized 
domain is denoted by The length e is proportional to the defect core size [21] and we assume 
that e <C The boundary, consists of four straight edges, Ci... C 4 , and four curved arcs of 

radius e enclosing the corners denoted by 71... 74 respectively. We use Green’s Theorem to write 
the energy as 


E = ^ [[ M\^dA = 



(63) 


where v is the outward pointing normal to each boundary segment. The contributions from the 
curved arcs, 7 *, are all of order e. We define four functions related to the integral contributions from 
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the straight edges; 


Si 


{N,b) :=8f2 


coth (^(271. — 1 ) In 6 ) + 1 


n=l 

oo 


2 n- 1 


, , csch (^( 2 n — 1 ) Inft) 

^ 2„._i ' 


coth (ynln 6 ) + 1 


n=l 

oo 


S3{N,b) :=8Y, 


n=l 

oo 


n 


s^iN,b) :=8j] 


csch (yn In 6 ) 


n=l 


n 


As e — 7 - 0 , the regularized energies within are then given by 


E r. KTr{log{^] +E] +0(e), 


(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 


where the normalized energy, E, is the interior distortion energy and the logarithmic term originates 
from the defects. The normalized energies of the four states are: 


Eui = 


si{N,b) + S 4 iN,b)-S 2 iN,b)-S 3 {N,b) , (iV + 2) 


+ 


41V 


1 


Eu 2 


si{N,b) + S 4 iN,b)-S 2 iN,b)-S 3 {N,b) , (iV-2) 


+ 


Eua = - 


siiN,b)+S2{N,b) , 1 


N' 


+ log r + n log T79 h 


41V 

166 




1 V 2 


Ed = 


S 2 {N,b)-si{N,b) , 1 


N 


+ — log - + - log ^ . 


166 




(69) 

(70) 

(71) 

(72) 


We make some immediate comments. In the case of a square or rectangular domain, the curved 
arcs, 7 j, make an order contribution to the OF energy M- In the case of an annulus, the curved 
arcs ji make an order e energy contribution due to the curvature of the boundary. 

In Figure we plot the normalized OF energy of the four different states, as a function of IV, in 
a fixed sector 11for some fixed values of 6. For small values of IV, the U 2 state is the minimum 
energy state in the set {Ui, U 2 , U 3 , D}. In contrast, for a square or rectangle, the diagonal state has 
the minimum normalized energy |14) . As N increases, there is a cross-over between the diagonal 
and rotated states and the D state has the minimum normalized energy state for large IV. This 
is consistent with the fact that Ow approaches a rectangle as IV —>■ 00 . The cross-over critical 
value, N = N^, depends on 6 ; in particular, our limited simulations suggest that Nc is an increasing 
function of 6 . This gives insight into preferred defect locations as a function of the annular aspect 
ratio. 

In all cases, the defect-free state has lower energy than the generalized diagonal or rotated 
solutions, {Ui, U 2 , U 3 , D}, for small values of e, see Figure]^ This stems from the A'log(l/e)- energy 
contribution from defects, which is absent in the defect-free case. 


5.2 Elastic Anisotropy 

We use the commercial PDF solver, COMSOL, to numerically solve the Euler-Lagrange equation 
Q with S ^ 0. We find all 4 states, {Ui, U 2 , Us,D}, within a sector see Figure]^ We further 
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Figure 7: The normalized energy of the four states plotted against N, for fixed values of b. 



Figure 8: The one constant Oseen-Frank energy of the three lowest energy states; the defect-free 
state and the U 2 states with = 1, 2, for defect core size e = 0.002. 

plot the structural differences between the one-constant case (with 5 = 0) and the anisotropic cases 
in Figure The structural details in the director are qualitatively the same. 

The regularized OF energy within can be expressed as the sum of the defect contributions 
and the normalized energy, E( 6 , b, N). 

i^~X3vr((l-^)logQ) +Ei6,b,N)^ +0(6). (73) 

In Figure [Toj we plot the numerically evaluated normalized energies for fixed b and N. For small N, 
the rotated U 2 state has the minimum energy in the set, {Ui, U 2 , U 3 , D}, for all admissible values of 5 
and b. As 6 —)• 1, the normalized energies oi Ui,U 2 , D converge whereas the normalized energy of the 
U 3 state diverges as 6 —)• 1. From Figure the Os-director is constrained to rotate by approximately 
vr radians over a distance of (1 — 6) units, leading to the energy blow-up as 6 —)■ 1. 

For suitable values of e and 6 close to 1, the defect-free state ceases to be the global minimizer of 
the Oseen-Frank energy, see Figure pT| As b decreases, we see the U 2 state with N = 1 minimizes the 
energy, than as b decreases further the U 2 state with N = 2 becomes the minimizing director field. 
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Ug state, 5 = -0.5 




Ug state, 8 = -2 




Figure 9: Examples of the effect of elastic anisotropy on minimizers of the Oseen-Frank energy in 
^2 and b = 0.25. We show numerically calculated director fields (left) and the difference in radians 
between the numerically evaluated 6 and the one constant solution given in Equation (58) (right). 


This cross-over in energies demonstrates the importance of elastic anisotropy and has implications 
for experimental observations, m. 

6 The Defect-Free State in the Landau-de Gennes Model 

This section focuses on the defect-free state in the two-dimensional LdG theory. We study maps 
Q : n —7- ^2 where = |xGM^:6<|x|<1} and S 2 is the space of symmetric traceless 2x2 
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6 = 0.1 
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6 = 0.3 



Figure 10: The normalized energies of the four states in Q 2 , for fixed values of b. 



Figure 11: The Oseen-Frank energy, with 5 = 0.9, of the three lowest energy states; the defect-free 
state and the U 2 states with = 1, 2, for defect core size e = 0.002. 


matrices, subject to the fixed boundary condition 


Qb 


-(n®n 


m (8) m) 


(74) 


on r = 6 and r = 1. In ( fTi] ), n = (cos (j), sin (/>), m = (— sin (j), cos 0) and cj) is the standard polar angle 
defined by taniji = We define the defect-free state to be 


Q* = s{r) (n (8) n — m (8) m), 


(75) 


where s : —)■ M is the unknown scalar order parameter subject to the Dirichlet conditions, s{b) = 
s(l) = The LdG energy of Q* can be computed using (3) and (4) to be 


I[Q*] = 27r 




(76) 


and we define the optimal order parameter s to be a minimizer of ( |76[ ) subject to the Dirichlet 
conditions on r = 6 and r = 1. It is straightforward to verify (see |25l |26] ) that the minimizing s is 
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a classical non-negative solution of the following second-order ordinary differential equation 


,"+■^-^ 4 = - 1 ). 

r L 


( 77 ) 


subject to s{b) = s(l) = If s were negative for tq < r < ri, then one could define the admissible 
comparison map, w = |s| which would have the same energy as s i.e. I[w] = /[s] but would 
have discontinuous derivatives at points where s vanishes, contradicting the regularity properties of 
minimizing solutions. One can check that Q*, thus defined, is indeed a critical point of the LdG 
energy in ([^ and the second variation of the LdG energy about Q* is given by 

:= jj ^|VV|2 + ^(2s 2 _ i)|v|2 + ^(Q* • VfdA, (78) 


where V(r, cj)) G S 2 is an admissible perturbation about Q* with V(6, cj)) = V(l, cp) = 0. 
We use the following basis for the space of symmetric, traceless 2x2 matrices, 


and 


E = 


F = 


Then, we can write an arbitrary V as 


sin(2(/)) \ 

— cos{2(f>) J 

(79) 

cos(2iji) \ 
sin(2(/)) J ' 

(80) 

-w{r, (/))¥, 

(81) 


where v{b,(p) = v{l,(p) = w{b,(p) = w{l,<p) = 0 and v,w are 27r-periodic in 4>. We use a basis 
decomposition for the functions, v and w, as shown below: 

00 

v(r,4,) = Y,a n{r) cos{n(p) + hn{r) sm{n(p), (82) 

n=0 
00 

w{r, 4>) = c„(r) cos{n(j)) -|- dn{r) sm{n4>), (83) 

where an = bn = Cn = dn = 0 on r = b, 1 for all n. The second variation can, thus, be expressed as 

00 

(5 I[Q, ] — 27tLq[(1{)^ CIQ; Cq, Cq] ~\~ 27r ^ ^ Lfiio^n^ bfif C 725 dn]? 

n=l 


where the functional Ln is defined to be 

Ln[a, 6, c, d] := j ^(a')^ (6')^ -h (c)^ -h (d')^ -h ^ + b‘^ + c^ + d^) 

{ad — cd)^ r + — (a^ -|- 6^ -|- -|- d^) (2s^ — l)r -|- /^2 _|_ j^ 2 \ 

JL ' L ' 

(85) 
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Proposition 2i_ The defect-free state, Q*, defined in (|75|), is a locally stable equilibrium of the 

fe2 + l 

~2fF' 


LdG energy in (3) with r(;(Q, VQ) = 4|VQp, for all 0 < 6 < 1 and ^ 


Proof. We follow the same strategy as in where the author establishes the local stability of 
the defect-free state on a disc, as opposed to an annulus in two dimensions. For a disc, the order 
parameter is defined for 0 < r < 1, the order parameter vanishes at r = 0 and is a monotonically 
increasing function subject to a Dirichlet condition at r = 1. The function s, as defined in (77), 
cannot be a monotone function for 6 < r < 1 by virtue of the imposed boundary conditions, < | 
by an immediate application of the maximum principle and hence s must have an intermediate 
local minimum m- The key ingredient is to show that minLn[an, bn, Cn,dn] > 0 for non-trivial 
an,bn,Cn, dn and min Lo[ao, UO) cq, cq] > 0 for non-trivial oq and cq. 

Lemma 1: The second variation, > 0 if minLo > 0, minLi > 0 and minL 2 > 0. 

Proof. It suffices to show that for re > 1, 


Ln+ 2 {a,b,c,d) - Ln{a,b,c,d) = j + b^ + c^ + d^) + ^{ad- bc)^ rdr > 0 


Using Young’s inequality. 


(4re -|- 4)(a^ -I- 6 ^ -|- -|- 16(ad — be) > (4re — 4)(o^ + b‘^ + + d)^) >0 


( 86 ) 

(87) 

□ 


The problem of stability now reduces to a study of the functionals Lq, Li and L 2 . The integral 
Lo is given by 


Ln 


^ +^{a‘l + 4)r dr +- l){al +cl)r+ ^s‘^alr dr. {88) 


It suffices to show that Lq > Iq = ^ + ^~T “ l)agrdr > 0 for all admissible oq such 

that ao{b) = ao(l) = 0. We define a{r) := s(r)ao(r) and recall the governing ordinary differential 


equation for s in (77). Then 


rdfds^ 

In = / u- s—a r s 


L dr \ dr 


rdr = 

dr J 


da , 

— I rdr > 0 
dr 


(89) 


for any non-trivial a : [ 6 , 1 ] —)■ M with a{b) = a(l) = 0 . 
It is straightforward to check that 


L^ 


1 (01,61, Cl, di) > Li( \Jal + 6f,0, ^Jc\ + d\,Q 


(90) 


by using the change of variable, oi = y^of -|- 6 f cos d, 61 = \/a\ + b‘\ sin 9 in the definition of Li 
above. Let a = af + bf and d = ^/~(^~+~d^ where a{b) = a(l) = d( 6 ) = d(l) = 0 . At this step, we 
diverge from We minimize the functional ( [76| ), subject to the boundary conditions s{b) = 0 and 

s(l) = and define re to be a corresponding energy minimizer. The existence of such a function re 
is guaranteed from the direct methods in the calculus of variations and the minimizer re is a classical 
solution of 

o/ All \ A\ 

( 91 ) 


,, re' 4re 2 

re -I-^ = —(2re - l)re 

r r^ L 


with re( 6 ) = 0 and re(l) = 
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Proposition 3: The function u is unique, monotonically increasing and u < s for all r. 

Proof. One can show that u is unique and monotonically increasing by an immediate adaptation 
of the arguments in |26j . which are omitted here for brevity. We show that u < s. We assume 
for a contradiction that s(r) < u{r) for some r G {b,l). Then, there must exist ri and r 2 , where 
6 < n < r 2 < 1, for which s(ri) — u{ri) = 0, s'(ri) — u'{ri) < 0, s(r 2 ) — u{r 2 ) = 0, s'(r 2 ) — u'{r 2 ) > 0 
and s(r) — u{r) < 0 for r G [ri, r 2 ]. We multiply the differential equations for s and u by ru and rs, 
subtract and integrate over (ri,r 2 ) to find 


rr 

J r^ 


rus" — rsu" + us' — su'dr = 


L 


us{ 2 s^ — 2 u^)rdr. 


(92) 


The RHS is negative since s — u < 0. The LHS is 

|•r2 r2 

LHS = / u{rs')' — s{ru')'dr = {rus' — rsu') 

J ri n 

= r 2 s{r 2 ){s'{r 2 ) — u {r 2 )) — ris(ri)(s'(ri) — u'(r 2 )) > 0, (93) 

yielding the desired contradiction. □ 

Since s > u, we immediately have 

Li(a, 0, d, 0) > J ^(a')^ + (d')^ + ^ (a^ + d^) - r 

+ ^ (a^ + d^) {2u^ — l)r + ^i^a^u^rdr = m. (94) 

ij 1 j 

It remains to show that minm > 0. The rest follows by analogy with HZ! with some technical 
differences. We compute the Euler-Lagrange equations associated with a and d: 

,, a' 5a 4d 1^1/2 

a H- ^ = ——aiou — 1) — ammm, 

r r^ r-^ L 

d' 5d 4a |H| ,,9 . , 

d 4-H— 7 ^ = —-ld( 2 a^ “ 1 ) “ dmmm, 

where minm acts as a Lagrange multiplier. This system is satisfied weakly, with minm = 0 and 
no boundary conditions, by the functions A = u' and D = 2u/r. Multiplying the Euler-Lagrange 
equations (95) by rA and rD respectively and the system for A and D by ra and rd respectively 
followed by integrating over (6,1) yields 


— minm / r{aA +dD)dr = {r{a!A +d!D)') 


(95) 


The RHS is negative since a, d > 0 , a'(l),d'(l) < 0 and a'{b),d'{b) > 0. Therefore, minm > 0. 
If minm = 0, then we must have (a'(l)H(l) + y/2d'{l)) = ba'{b)A{b). This implies d'(l) = 0 and 
a'{b)u'{b) = a'(l)a'(l) = 0. For a'(l)u'(l) = 0 to hold, we must either have a'(l) = 0 (and hence 
a = d = 0) or u'{l) = 0. From the ODE for u in (91), we see that u"{l) = 2\/2, which implies 
that n(l) is a local minimum contradicting the fact that u is a monotonically increasing function for 
6 < r < 1. 


The last step is to show that L 2 is positive for large values of We first show that L 2 > 0 i 
00 by following the strategy laid out in |28) , and then derive an explicit condition 


m 


the limit as ^ 


in terms of b and We note that 


L2 > 


8(a2 - d2)" 


a '2 + d '2 + 2 ^ “ 1) + d\{ 2 s^ — 1))^ r dr = mi 


(96) 
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where and are positive functions of r. 

The first stage is to show that for sufficiently large mins(r) is controlled by 


— s ■ (2s'^ ■ - 11 + > 0 

L ri 


Rearranging, 


Proposition 4: ^ ^ ^ uniformly as ^ oo. 

Proof. The function s attains its minimum at r* G (b,!)- We define s(r*) = Smin- Therefore, 
s'(r*) = 0 , s"(r*) > 0 and 

1/11 /I/. . 

(97) 

(98) 

□ 

(99) 


/ N /l 2L 

S{r) > Smin > \ - - 


1 |R| 

14 10 ^ ^ —y oc 

2 |T|r2 V2 L 


We let a 2 (r) = s(r)a(r), d 2 {r) = s{r)d{r), and use Hardy’s trick to write mi as 


2 f 12 j/2 + <7^) ISad 4|H| 

mi= i a'2 + ^5- - -^ ^ 

. h V L 


and consider the terms Then 

4(a^ + d^) 16ad\ 4|H| 4 2 ^ + d?) 256a^ 4|H| 4 , 

^2 r"^ 1 ^ L L 1 

The coefficients of are positive if the order parameter Smin > | for all 6 G (0, !)• From Proposition 
4, we can always find some Ti such that Smin > | for all ^ > Ti. 

The coefficients of are positive if ^ > F 2 , where we dehne T 2 = ■ Then L 2 > 0 if 

^min^ 

> max{Ti,r2}. 

Finally, we show that L 2 is positive if ^ is greater than an explicit positive quantity dependent 
on b. We let a(r) = a(r) cos(0(r)) and d{r) = a(r) sin(0(r)); then 


, 2 1/2 8a^sin(20) 4|H| 22 2 /n\ \ , 

mi > I a ^ H- 2 - 2 - - -sim(0) ] r dr = m 2 . 


L 


( 101 ) 


The minimum of m 2 , with respect to 9, occurs at 


1 

= - arctan 
2 


4L 


\A\r‘^s{r^‘^ 


( 102 ) 


and 


/ 2 /2 I 

m 2 = s a + 


4a^ 

J.2 


32|H|a2 


r^Ls^iJl + 

We note that a = 0 is always a critical point of 

32|H|a2 


+ 2 ,V^ I 1 

Jj 


4a^ 

^2 


r^Ls^yJl + 


16L2 
I Aps^r^ 


+ 2 sV^ I 1 - 

IJ 


1 + 




1 + 


16L2 

|j4Ps4^4 


r dr. (103) 


(104) 
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(105) 


and is a minima if ^ > 


3 

^2 c2 ■ 


Using the fact that s > (see [2H]), then L 2 > 0 if 


1^1 . 3 ( 6 ^+ 1 )^ 

L 264 


□ 


7 Conclusions 

In this paper, we investigate nematic equilibria confined to two-dimensional annuli with strong or 
weak tangential anchoring modelled in the continuum Oseen-Frank and Landau-de Gennes theories. 

We compute explicit stability criteria for the ddefect-free state with a radially-invariant direc¬ 
tor field in both the Oseen-Frank and Landau-de Gennes theories. In the Oseen-Frank theory with 
strong anchoring, we demonstrate that the defect-free state undergoes a supercritical pitchfork bifur¬ 
cation and compute a new energetically preferable defect-free state with spiral-like director profile. 
Furthermore, we demonstrate that this spiral-like state is locally stable for <5 —)• 1. In the Oseen- 
Frank theory with weak anchoring, we produce a new stability diagram (Figure in terms of the 
parameters 6 and a. We demonstrate that the defect-free state is stable with respect to symmetry¬ 
breaking perturbations when the outer radius exceeds a critical value, which is precisely the surface 
extrapolation length of the nematic material in use. In parallel, we obtain an explicit local stabil¬ 
ity criterion for the defect-free state, in terms of the temperature and geometry, in the Landau-de 
Gennes framework. 

We model nematic equilibria with defects in the Oseen-Frank theory on the boundary by com¬ 
puting local solutions of the Laplace equation on an annular sector with Dirichlet tangent boundary 
conditions. We find three rotated and one diagonal solution, by analogy with similar work done 
for squares and rectangles in [laini [30]. We compute analytic expressions for the corresponding 
director fields and their one-constant Oseen-Frank energy and these energy expressions allow us to 
distinguish between the energy contributions from the defective corners and the bulk distortion and 
give information about the optimal number and arrangement of boundary defects as a function of 
geometry. In particular, we find that the rotated JJi state, which has been found experimentally in 
mE], has the minimum energy in the restricted class {Ui, 1 / 2 , t/ 3 , 11 }, for small numbers of defects. 
In contrast, the diagonal solutions always have minimum energy for a square or a rectangle in the 
one-constant framework, see [T3|. Furthermore, we investigate the effect of elastic anisotropy and 
demonstrate that for iLs ^ Ki, the U 2 solution with boundary defects can have lower energy than 
the defect-free state for specific choices of the geometry. The singular limit <5 —>■ 1 will be investigated 
in future work. 
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